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Abstract 

A fit of amplitudes to the experimental branching ratios to two mesons 
is used to construct a new estimate of neutral D mixing which includes 
SU(3) breaking. The result is dominated by the experimental uncertain- 
ties. This suggests that the charm sector may not be as sensitive to new 
physics as previously thought and that long-distance calculations may not 
be useful. 
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Introduction 

The prospects of probing new physics in D°-D° mixing has been a topic recent 
discussion |3J p|. In this letter we will present a new estimate of Anic due 
to mixing via two-body hadronic intermediate states. Because the couplings 
between the D mesons and the possible intermediate states are not all known, 
we will rely on the results of a fit to the branching ratios ||. This fit utilizes 
a complete parameterization of the decays in the framework of broken flavor 
SU(3) and allows us to extract the (complex) couplings needed to estimate 
Amo. Intermediate states from the fit are limited to pseudoscalar-pseudoscalar 
(PP), pseudoscalar- vector (PV), and vector- vector (VV). We will find that the 
results are overwhelmed by their uncertainties. These uncertainties owe their 
origin to the large (sometimes ~30%) uncertainties on the experimentally mea- 
sured branching fractions. Other estimates, based on short-distance processes 
or on the heavy-quark effective theory, give Am^ to be two orders of magnitude 
smaller that our central values and uncertainties. In the absence of cancellations 
among the hadronic modes, our estimate allows the standard-model contribu- 
tions to D°-D mixing to be close to the current experimental limit ~ 10 



1 Formalism 

Weak interactions to second order in the coupling Gf give an off-diagonal part 
to the Hamiltonian which represents D° <-» D transitions. Such a term is 
responsible for the mixing. Its real part enters the mass matrix and generates 
a mass difference between the two eigenstates. In the absence of CP violation, 
we can write this as 



GeVl. 




(1) 



The mass eigenstates are 




(2) 



with masses 



mi = M + |Am, 
m 2 = M - |Ato. 



(3) 
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The mass difference Am is a convenient parameter for the size of the mixing. It 
is the quantity whose measurement and calculation are important in the question 
of probing new physics in D°-D° mixing. 

Before continuing, we must draw the distinction between short- and long- 
distance estimates to Am. Short-distance contributions come from the calcu- 
lation of the "box" diagrams in the quark model. Long-distance contributions 
come from the dispersion due to intermediate hadronic states. In this case, the 
approach is only valid for internal momenta that are below the scale at which 
QCD becomes nonperturbative. Typically this scale is taken to be 

fi ~ 1 GeV. (4) 

The dependence on \i in our estimate will be explicit. 



2 Modelling the Couplings 

In order to construct our estimate, we need to know the couplings D°MiM 2 
between the neutral D mesons and the lowest-lying nonets of pseudoscalars 
and vectors. To this end, we have parameterized these couplings in the SU(3) 
framework with flavor-symmetry breaking by an octet. The details can be found 
in 0. Here we will only sketch the calculation. 

The particles are organized into SU(3) multiplets. The D mesons (D°, D + , 
D s ) form an antitriplet 3. The pseudoscalars (n^, 7r°, K , K°, K , T], t/) fit 
into an octet and a singlet. The physical r\ and rf are mixtures of the octet and 
singlet pieces, with mixing angle given by experiment 0. The vectors (p , p°, 
K*^, K*°, K*° , u, (j)) are likewise arranged, with a vector mixing angle that is 
found in [|J. 

The charm-changing Hamiltonian is proportional to a product of quark 
currents 

Hac=i ~ - 75)<? - 75)c, (5) 

where 

q = cos Bed — sin Oqs, 
q' = sin 9 c d + cos 6 c s. 

Since the quark-annihilation operators in H transform as antitriplets, and the q 
operators as triplets, we can expand H in terms of SU(3) representations. The 
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Hamiltonian is found to transform as 15 and 6. The Clebsch-Gordan factors in 
this expansion and in the following were calculated by computer || 0. 

The amplitude for each decay of the type D — > PP, PV, VV can now be 
written as the sum of reduced matrix elements with appropriate Clebsch factors. 
For each of PP, PV, and VV there are 48 (complex) reduced matrix elements. 
The number of parameters is too large to be fit by the available 45 measured 
modes and 13 modes with experimental limits ||. Therefore, we must make 
some assumptions to limit their number, first, we assume that corresponding 
reduced matrix elements of PP, PV, and VV are proportional in magnitude. 
This proportionality is represented by two new parameters A PV /pp and A/v/pp- 
Second, we assume that the phase of each reduced matrix element is given by 
the phase of the representation into which the daughter mesons are contracted. 
These are the {mvi)i, {ViP)b, (PP)i, (PP)s, (PP)a7, (Wi)i, faiV) 8 , (^P) 8 , 

(pv)!, (pv) 8 , (pvv, (pv) 10 , (pv)io, (pv) 27 , (u^Oi, KV) 8 , (w)i, (vv) 8 , 

(VV) 2 7- The phases then become new parameters. With these assumptions, the 
number of linear combinations of reduced matrix elements that contribute to 
any decay is reduced to 40. 

There are now far fewer parameters than we began with. The data constrain 
all of them, except for three combinations of reduced matrix elements and the 
phases of (771/71)1, (^1^1)1, and (771^1)1. There are too many free parameters in 
the singlet-singlet cases, and so we do not attempt to make any estimates of 
their values. However, one of the remaining combinations of matrix elements 
that involve the D° can be constrained by reasonable estimates of one additional 
decay mode. In order to see the effect of our lack of knowledge in this case, we 
make two different estimates, called schemes A and B. In the former we use 

B(D° — > r]K°) = 3tan 4 # c B(D° — > rjK°), (7) 

and in the latter we take 

B(D° — > 4>K°) = 3tan 4 # c B(D° — > (fK°). (8) 

The coefficient 3 is motivated by the size of the recently measured mode D° — > 
K + 7v~ 0. Interestingly, a coefficient of 1 does not result in a consistent fit. 
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3 Estimate of Am 



The long-distance contributions to Am arise from dispersive effects involving 
intermediate states. The two-body hadronic intermediate states were considered 
by [IIJ are were estimated by 

Am K± ^ = ^ln^f [T(D° ^ K+K-) + T(D° ^ n+iT-) 

-2^T(D° ->• K+7i')T(D° -»• AT-7r+)l . 

Here [i is the cutoff discussed in Section 2. Notice the implicit assumption that 
the couplings are relatively real. If we insert the most recent values for these 
rates H ||, we obtain 

Am K± ' n± = (-0.75 to 0.29) x 10~ 15 GeV. (10) 

The range of values is due to the uncertainty on the branching fraction to K + ty~ 
||. The purpose of this exercise is to show that large SU(3) breaking may 
give large long-distance contributions to neutral D mixing. Should the SU(3) 
breaking be due only to the K-n mass difference, the rates would be related by 

T(D° -> K + K-)$(KK) = T(D° -► 7r+7r-)$(7T7r) 

= tan 2 6 c T(D -> K-n+)$(Kn) (11) 
= cot 2 e c T(D° -> K+n-)$(Kn), 

where $ represents the phase-space corrections. These corrections are on the 
order of a percent, and so in this case 

Am ~ 10" 17 GeV, (12) 

a value consistent with other estimates, as discussed below. 

We will make two improvements on this approach. First, we will include all 
PP, PV, and VV intermediate states, with the exception of the singlet-singlet 
states. Second, we will allow the couplings to take complex values. Equation |9| 
is replaced by 

Am L D D = — In ^ x A(D° -> I)A*(D° -> 7), (13) 

where N is a normalization factor given by 

T(D° ^ K-ii + ) = N\A(D° ^ K-n + )\ 2 . (14) 



4 



The coupling are extracted from the fit of the previous section. 

Our estimates for Am are presented in the Table. The estimates due to PV 
intermediate states are consistent with zero and have uncertainties on the order 
of 50 xlO -15 GeV. The estimates due to PP and VV intermediate states vary 
according to our choice of estimate for the unconstrained modes (Equation |7| or 
H) and are at most one standard deviation from zero. The source of the large 
uncertainties is the uncertainty with which we know the individual branching 
fractions used to constrain our parameterization. The contributions from indi- 
vidual modes enter with differing phases, but the uncertainties are cumulative. 
The result is that our attempt to estimate Am is nearly overwhelmed by un- 
certainties. It is also worth noting that these uncertainties are all greater than 
the prior estimate which used only the modes n + n^ , K + n~ , K~n + , and K + K~ 
(see Equation |ToD . 



4 Discussion 



We have found that a long-distance calculation of Amp is dominated by the 
experimental uncertainties. The treatment of D-meson mixing by Donoghue et 
al. JT(| has underestimated these uncertainties. On the other hand, there also 
exist estimates based on the underlying quark processes. The short-distance 



calculation based on the "box" diagrams has been done in JTT] [12 [13] [10| |L4 



For example, |TTJ] finds that 
G F 



Am 



box 



O 



v/2 4tt sin 2 6 v 



cos 2 9c sin 2 8, 



m; 



ni 



2\2 



c~ 



M^ml 



-m D F 2 D {B D -2B' D ). (15) 



Here B and B' are hadronic factors defined in |L0[ and Fd is the pseudoscalar 
decay constant: 

(0\A^\D°(p)) = iF DPfl . (16) 

If we assume that the hadronic factors B ~ B' ~ 1, and take Fd ~ 300 MeV 
0] and m s ~ 250 MeV, then the contribution of the box diagrams to Am is 



Am 



box 



4.8 x 10~ 17 GeV. 



(17) 



In addition, an estimate of the mass difference in the heavy-quark effective field 



theory (HQEFT) was performed in [15 | 16|| . These arrive at estimates for the 
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contributions of 4-, 6-, and 8-quark operators: 




1 xl(T 17 GeV, 

2 xlO~ 17 GeV, 
0.5 xl(r 17 GeV. 



(18) 



These estimates replace both the long-distance calculation and the box calcu- 
lation. In order to reconcile the HQEFT estimates with our approach, there 
must be cancellation of the individual long-distance contributions to about one 
percent. In addition, it may be that the HQEFT estimate has underestimated 
the SU (3) breaking involved in the D system. The breaking of SU(3) is known 
from the splittings of hadronic masses to be at the level of 20-30%. However, in 
the charm system, we know from the branching fractions that SU (3) is broken 
at the 100% level. In the HQEFT estimate, only the quark-mass differences 
were used to break the flavor symmetry. We conjecture that this is the reason 
that the HQEFT estimate is on the same order of magnitude as the box and 
much smaller than the long-distance estimates. 

The cancellation between contributions from different hadronic interme- 
diate states needed to reconcile the HQEFT approach with that of the long- 
distance estimates would have to be among, rather than within, the individual 
SU(3) representations. We can see this by considering the contribution due to 
the complete octet of pseudoscalar mesons. Although the singlet-singlet pieces 
of the PP decay modes are completely unconstrained and the singlet-octet only 
partially constrained, the octet-octet parts of the amplitudes are completely de- 
termined by data. Therefore we are able to extract the octet-octet parts from 
the fit without relying on the assumption of Equation [7] or |8|. It is then possi- 
ble to construct an estimate of Am that includes the entire pseudoscalar octet, 
without the mixing in of singlet pieces. We find that this contribution is 



This estimate differs significantly from zero, and is therefore an indication that 
cancellations among hadronic modes must be between the various SU(3) repre- 
sentations. 

It is not possible for us to determine whether cancellation occurs within 
the complete set of PP intermediate states, including both SU(3) octet and 
singlet pieces. The expected size of the PP singlet-octet contribution can be 



Am 



P octet 
D 



(9.6 ±2.2) x 10~ 15 GeV. 



(19) 
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seen in the difference between Equation [TP] and the entries in the first column of 



Table 1. These singlet-octet contributions vary by our estimation schemes for 
the unconstrained modes, and are on the order of half of Am^> octet . We have no 
expectations on the size of the singlet-singlet contribution, but can remark that 
it would be entirely due to SU(3) breaking. It is possible that PP modes will 
cancel among themselves, once the singlet- singlet and singlet-octet pieces are 
fully included. But should that not occur, then the only hope of cancellation 
would be between the PP, PV, and VV modes taken together. 

Due to the large uncertainties (sometimes ~ 30%) on the experimentally 
measured branching fractions, the uncertainties on our estimates of the long- 
distance contribution are very large. There are two things that could be learned 
from this. First, the large uncertainties indicate that we are yet unable to 
make a useful calculation of Amp using intermediate hadronic states. It is 
not reasonable for the experimental situation to improve to the point where the 
uncertainties on Am^ reach the precision of the box or HQEFT estimates in the 
near future. Continued endeavors using one-particle (resonant) or three-particle 
intermediate states are discouraged. Second, in the absence of cancellations 
among the hadronic intermediate states, the large range of possible values allows 
the standard-physics contributions to D-meson mixing to be as large as 

\Am D \ = l(T 13 GeV. (20) 

This is near the experimental limit [|J, and therefore this process would be less 
likely to be useful as a probe of new physics than previously thought []J . We 
are left with the hope that a direct measurement of the D-meson mass difference 
can be obtained and that these questions can be resolved. 
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Table : Estimates of Amp. PP, PV, and VV refer to pseudoscalar- 

pseudoscalar, pseudoscalar-vector, and vector-vector intermedi- 
ate states. Schemes A and B for estimating some doubly sup- 
pressed neutral modes are explained in Section 2. The first line 
neglects those modes. The last line includes the full pseudoscalar 
octet. In it, the octet-octet parts of 7777, r]7]', and rfrf are included 
and the singlet-octet parts are excluded. All values are in 1CT 15 
GeV. 





PP 


VV 


PV 




no estimates 


7.3 ± 3.9 


19.1 ± 11.3 


-60.3 ± 


63.3 


scheme A 


10.1 ± 4.4 


25.5 ± 11.9 


-56.5 ± 


63.9 


scheme B 


4.6 ± 4.5 


14.7 ± 12.1 


-65.5 ± 


63.9 


K ± and 7r ± 


3.7 ± 1.3 








full octet 


9.6 ± 2.2 
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